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ABSTRACT 



Aims. The aim of this paper is to provide a measurement of the mass accretion rate in a large, complete sample of objects in 
the core of the star forming region p Oph. 

Methods. The sample includes most of the objects (104 out of 111) with evidence of a circumstellar disk from mid-infrared 
photometry; it covers a stellar mass range from about 0.03 to 3 M0 and it is complete to a limiting mass of ^^0.05 Mq. We 
used J and K-band spectra to derive the mass accretion rate of each object from the intensity of the hydrogen recombination 
lines, Pa/3 or Br7. For comparison, we also obtained similar spectra of 35 diskless objects. 

Results. The results show that emission in these lines is only seen in stars with disks, and can be used as an indicator of 
accretion. However, the converse does not hold, as about 50% of our disk objects do not have detectable line emission. The 
measured accretion rates show a strong correlation with the mass of the central object {MaccOC Ml'^^°'^) and a large spread, 
of two orders of magnitude at least, for any interval of M*. A comparison with existing data for Taurus shows that the objects 
in the two regions have similar behaviour, at least for objects more massive than ^0.1 Mq. The implications of these results 
are briefly discussed. 
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, ^ 1. Introduction 



Accretion disks are common around young stars of all 
mass, from M^.'^S Mq down to very low mass objects and 
brown dwarfs. They form during the collapse of the molec- 
ular core, from which the star is born, and last well beyond 
this initial phase, when the core has dispersed and the star 
has acquired most of its final mass. The accretion disks are 
the birthplace of planets, whose formation and evolution 
are controlled by the disk physics. 

Even if accretion disks have been part of the accepted 
paradigm of star formation for many years, many of their 
physical properties are poorly known, and the physical 
mechanism of angular momentum transfer, which deter- 
mines the disk evolution, is still unclear. The physical 
quantity that controls the accretion phase is the mass 
accretion rate through the disk Mace- This quantity can 
be derived only indirectly, by fitting models to observed 
quantities such as the UV excess emission and/or the pro- 
files and intensity of lines believed to form in the accret- 
ing gas. Measurements of Mace are now available for a 
large number of stars in Taurus (e.g., MuzeroUe et al. 120051 
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* Based on observations collected at the European Southern 
Observatory, Chile. Program 073.C-0179. 



and references therein). The results have shown that Mace 
is a strong function of the mass of the central object, 
roughly oc M^, and that a large dispersion is present 
(about two orders of magnitude) for objects with the same 
M^. Both results are a challenge for accretion disk models, 
as discussed, e.g., by MuzeroUe et al. H2003fl and Natta et 
al. PTOjl . 

Measurements of accretion rates in other star forming 
regions are scarce in comparison, mostly limited to very 
low mass objects (MuzeroUe et al. I20U3I I2005|l . In a study 
of very low mass objects and brown dwarfs in Ophiuchus, 
Natta et al. (2004) found that they are actively accret- 
ing with Mace higher by at least one order of magnitude 
than objects of similar mass in Taurus. This could be due 
to a difference in age, since the Ophiuchus BDs are very 
young objects, younger than their Taurus counterparts, 
but could also be due to different environmental condi- 
tions. 

While it is clearly necessary to improve the physical 
models of accretion disks, at the same time it is important 
to study large and if possible complete samples of stars in 
a variety of star forming regions, differing in age and global 
properties. 

We report in this paper the results of a project aimed 
at measuring the mass accretion rate of a large sample 
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of pre-main sequence objects, ranging from a few solar 
masses to few tens of Jupiter masses, in the star form- 
ing region p Oph. The core of p Oph is perfectly suited 
for such a study, as it is rich in pre-main sequence stars, 
which include intermediate mass objects, T Tauri stars 
(TTS) and brown dwarfs (BDs). Its stellar content has 
been studied, e.g., by Luhman & Rieke (11999': LR99 in the 
following, and references therein to previous work), Natta 
et al. H2Ul)2fl and, more recently, by Wilking et al. H2UU5|I . 
Moreover, p Oph is very different from Taurus, younger 
and more compact, and it will allow us to explore the ac- 
cretion properties of pre-main sequence stars under differ- 
ent conditions, following the results of Natta et al. H2()()4L 
Ophiuchus has been observed in two mid-IR bands 
with ISO by Bontemps et al. 1^711711 BKAOl in the fol- 
lowing), who detected 199 sources in the p Oph core. Of 
these. 111 were classified, on the basis of their IR colors, 
as Class II objects, i.e., visible young stellar objects with 
evidence of disks. They provide a sample of systems with 
disks complete to a limiting mass of about 0.05 M0. In 
a spectroscopic study of the very low luminosity objects 
of the BKAOl sample, Natta et al. H2()02|l confirmed that 
they were BDs with mid-IR excess, very likely from a cir- 
cumstellar disk; as mentioned, these BDs show significant 
differences in accretion properties from their analogs in 
Taurus. 

The disadvantage of observing p Oph is its high extinc- 
tion, which makes veiling measurements in the UV and vi- 
sual impossible except for a few objects. The most effective 
way to determine Mace for the Ophiuchus sample is there- 
fore to use the luminosity of hydrogen recombination lines, 
such as Pa/3 and/or Br7. The relation between IR line 
luminosity and accretion luminosity, independently mea- 
sured from the UV excess, was established by Muzerolle 
et al. ((TF98) for TTS, and by Calvet e t al. for in- 

termediate mass objects. Natta et al. (|2004|l extended it 
to very low mass objects, where Mace was determined by 
fitting the observed Ha profiles with the predictions of 
magnetospheric accretion models. 

In this paper, we present the results of a spectroscopic 
IR survey of Ophiuchus objects. In Sec. 2, we describe the 
properties of the observed sample, which includes almost 
all (104 out of 111) the Class II objects and a subset (35 
objects out of 77) of the diskless systems (Class III), also 
from the BKAOl survey, that we will use for comparison. 
The observations, data reduction and method of analysis 
are discussed in Sec. 3. The results are presented in Sec. 4 
and discussed in Sec. 5. Sec. 6 summarizes our conclusions. 

2. Characteristics of the observed sample 

2.1. The sample 

The most complete survey of young stellar objects in the 
p Oph Main Cloud (L 1688) is that obtained in two mid-IR 
bands (6.7 and 15.3 ^m) with ISOCAM (BAKOl). Based 
on the near and mid-IR colors, the objects were divided 
in Class I (accreting protostars). Class II and tentative 



Class II (pre-main sequence stars with IR excess typical 
of disks, like classical T Tauri stars or CTTS), and Class 
Ill/tentative Class III (objects with colors typical of stel- 
lar photospheres, like weak- line T Tauri stars or WTTS). 
BAKOl estimate that their Class II sample of 111 objects 
is complete to a limiting luminosity Lj,~0.03 L©, corre- 
sponding approximately to 0.05 M©. The Class III sam- 
ple is only complete to ^0.2 L© (about 0.15 Mq). Note 
that not all the Class III objects have been confirmed as 
p Oph members. Barsony et al. H2005|l have recently con- 
firmed the accuracy of the ISOCAM results with ground- 
based 10/im observations of a large subset of the BKAOl 
sources. 

Our sample includes 104 of the 111 Class Il/tentative 
Class II objects (Class II for simplicity in the following) 
listed by BAKOl in the p Oph core. Most of the spectra 
(96) were obtained in the J band; the remaining 8, of ob- 
jects too weak in J, in the K band; one object has been 
observed at both wavelengths. As a comparison sample, 
we observed 35 of the 77 Class III and tentative Class III 
(in the following. Class III) objects, 31 in the J band and 4 
in K. The objects and their properties are listed in Tables 
C.l and C.2 1. 

2.2. Stellar parameters 

The stellar properties (i.e., spectral type, luminosity, mass 
and radius) of the BAKOl sample are well known only for 
a handful of objects. The main difficulty comes from the 
large uncertainties that affect spectral types, due to the 
combination of high extinction and large veiling, even at 
near-IR wavelengths (e.g., LR99 and references therein; 
Doppman et al. I^HMI Wilking et al. I^nn5|l . LR99, using 
K-band low resolution spectra, provide spectral types for 
37 of our Class II objects. However, 23 of them have un- 
certainties of almost one spectral class. 

Given the uncertainties, and considering that most of 
our objects do not have any spectral classification, we have 
decided to adopt a statistical approach, following BAKOl. 

First, we compute the extinction toward each object 
from the observed (J-H)-(H-K) colors, as given by 2MASS 
^, corrected to CIT system, adopting the Ophiuchus ex- 
tinction law of Kenyon et al. H1998|l and the locus of CTTS 
defined by Meyer et al. l|1997|l . The result can be expressed 
as: 

Aj = 2.31 (1.72 (J - H)ciT -{H- K)cit - 0.896) (1) 

This relation gives the correct reddening also for objects 
with no excess in the near infared, as long as they have 
Teff<5000 K (spectral type later than K2), which apphes 

^ All the tables are available in electronic form only. 

^ This publication makes use of data products from the 
Two Micron All Sky Survey, which is a joint project of the 
University of Massachusetts and the Infrared Processing and 
Analysis Center/California Institute of Technology, funded by 
the National Aeronautics and Space Administration and the 
National Science Foundation. 
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to the majority of our sample. However, for diskless ob- 
jects of earlier spectral type it will underestimate Aj and 
thus the inferred L^, significantly. We will come back to 
this point in §4.1. 

The stellar luminosity is computed from the J magni- 
tude and Aj, using a bolometric correction similar to that 
adopted by BAKOl: 

logi* = 1.24 + l.llogLj (2) 

where L^,, Lj are in units if Lq and Lj = 301 x 
-l^Q-(j-Aj)/2.5_ Eq,(2) assumes that the J-band disk emis- 
sion is negligible in comparison to the photospheric one. 
If this is not the case, the equation overestimates L^,. A 
recent work by Cieza et al. (|2UU5|I suggests that classical 
TTS have a J-band excess of ~ 0.3 mag on average, so 
that derived from the J magnitude is higher than the 
true one by about 30%. We have checked that, if such 
a correction applies to objects of all masses, none of our 
conclusions will change. Our values of L^, are in agree- 
ment with LR99 estimates always within a factor of 2, 
and generally much better. We have also compared the 
luminosities derived in this way with the results obtained 
by Natta et al. H2()()2|l from near-infrared J,H,K low resolu- 
tion spectroscopy for a group of 10 BDs; also in this case, 
the results are within a factor of 2, with the exception of 
one object (pOph-ISO 033), for which we underestimate 
the luminosity by almost one order of magnitude. This 
discrepancy has no impact on the results of this paper. 

When only H and K magnitudes were available (17 
Class II and 1 Class III objects), we estimated the stellar 
luminosity using eq.(2) and (4) of BKAOl. 

There are 5 objects (4 Class II and 1 Class III) that 
have companions clearly seen in our spectra, but which are 
not resolved in the 2MASS photometry. All the compan- 
ions have a good detection of the continuum; the flux ratio 
between the primary and the secondary is always larger 
than a factor of 3. Two of the companions (pOph-ISO 068b 
and pOph-ISO 072b) have been detected in the K-band 
p Oph multiplicity survey of Ratzka et al. H2005|l . with 
flux ratios to the primary of 0.19 and 0.16, respectively. 
We have accordingly not corrected the 2MASS magnitudes 
of the primaries for the contribution of the companions, 
because the corrections to the derived parameters would 
have been within the uncertainties. The secondary compo- 
nents have no detectable Pa/3 emission, and we will omit 
them from our analysis in the following; their properties 
are summarized in Table C.3. 

To determine stellar radii and masses, we make the 
assumption that the star formation in Ophiuchus is co- 
eval, and that all the objects lie on a single isochrone in 
the HR diagram. With this assumption, we can derive 
stellar mass, temperature and radius from the measured 
L*. This procedure is reasonable for the Ophiuchus core, 
whose age estimates range between 0.5 and 1 Myr, with 
very few stars older than that (BAKOl; LR99). In the fol- 
lowing, we adopt the D'Antona and Mazzitelli (,1997., and 
1998 web updates; DM98 in the following ) evolutionary 



tracks for an age of 0.5My. The uncertainties introduced 
by the assumption of coeval star formation and the differ- 
ences expected if other evolutionary tracks were used are 
discussed in Appendix A. 

The values of the stellar parameters are given in Tables 
C.l and C.2. 

3. Observations and data analysis 

3.1. Observations and data reduction 

Near infrared moderate resolution J and K band spec- 
troscopic observations of all targets in our sample were 
obtained at the ESO Observatories in Chile. The objects 
were either observed using the SofI instrument at the NTT 
3.6m telescope (June 2004, Visitor Mode) or the ISAAC 
instrument at the Antu 8.2m VLT unit telescope (Spring 
2004, Service Mode), as specified in Tables C.l and C.2. 
Detailed descriptions of both these instruments are avail- 
able on the ESO web pages ^. For all the objects that 
were observable at J-band, with SofI we used the 0.6 arc- 
second slit and the Blue low resolution grims, resulting 
in a spectral resolution of approximately A/AA ~ 1000 
and a spectral coverage from ^ 0.95 to ~ 1.64 /im; with 
ISAAC we employed the short-wavelength low resolution 
spectral mode with central wavelength 1.25 /im and 0.6" 
slit width, giving a spectral resolution of A/AA ~ 900 and 
a spectral coverage limited to the J-band. A number of 
objects were only observable at K-band, for these we ei- 
ther used the SofI Red low resolution grism with similar 
spectral resolution as for the Blue grism observations and 
spectral coverage from ^ 1.6 to 2.5 /im, or the ISAAC 
short wavelength low resolution mode with central wave- 
length 2.2 /im, which offers a similar spectral resolution as 
the J-band observations and a spectral coverage limited to 
the K-band. Integration times varied from about 0.5 to 2 
hours on source, depending on the expected brightness of 
the objects and observing conditions (in Visitor Mode). 

During the Visitor Mode observations at the NTT tele- 
scope, we acquired several telluric standard stars per night 
at varying airmasses; each Observing Block from our pro- 
gramme executed in Service Mode at the VLT was pre- 
ceded or followed by a telluric standard observed with the 
same instrument mode and at a similar airmass as our 
target stars. Spectroscopic flat fields and arcs were ob- 
tained during daytime either before or after our observa- 
tions. Standard methods were employed to calibrate our 
data. We did not attempt to obtain flux calibrated spec- 
tra; all our spectra are wavelength calibrated using OH 
airglow lines and corrected on an arbitrary intensity scale 
for telluric absorption and instrument response using the 
telluric standard star observations. 

Correction for telluric absorption and instrumental re- 
sponse was obtained observing at similar airmasses early 
type stars (early B or O) of known spectral type from the 

http:/ /www. Is. eso.org/lasiIla/sciops/ntt/sofi/and 
_http: / /www. eso.org/instruments / isaac / ^ 
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telluric standards lists of ISAAC ^. These stars all have 
Pa/3 or Br7 absorption which were manually removed from 
the spectra before applying the correction. 

Most of the spectra are of excellent quality; the detec- 
tion limits of the Pa/3 or Br7 equivalent width are in gen- 
eral of yhe order of 0.5-lA. Variations around this limit 
are mainly related to the signal to noise ratio achieved on 
the photospheric continuum of the individual objects. The 
signal to noise ratio depends on the telescope/instrument 
used, the observing conditions, the integration time and 
the apparent magnitude of the object. It is not necessar- 
ily a function of the object intrinsic luminosity because 
the extinction can be very different and because we tried 
as much as possible to observe two objects at the same 
time by properly aligning the slit, so that some relatively 
bright source near a faint one may have been observed 
with ISAAC and a long integration time. However, most 
of the lower luminosity objects have been observed with 
ISAAC and, expecting lower line intensities, with a higher 
signal to noise ratio; thus, the line detection limits for low 
luminosity objects are generally lower than for intermedi- 
ate luminosity ones. 

The sample studied in this paper includes also the 9 
BKAOl sources for which Natta et al. H2()()4|l obtained J 
and K band spectra with ISAAC. We have taken the Natta 
et al. l|^04l) J band spectra and reanalyzed them in the 
same manner used for the others. 

3.2. Method 

The luminosity of Pa/3 and Br7 are computed from the 
measured equivalent widths of the emission lines and the 
broad-band J and K fluxes, corrected for extinction, deter- 
mined as described in §2.2. No correction for underlying 
photospheric absorption was applied, since the expected 
equivalent width is small (<0.5 A; Wallace et al. |2000) 
for objects with Teff<5000 K, which represent the quasi- 
totality of our sample (see Table C.l and would not change 
the results. 

There are 12 Class II (11 of them have no Pa/3 detec- 
tion) for which it was not possible to determine line fluxes, 
due to lack of J magnitudes; they will not be included in 
the following discussion. Similarly, we will not consider 
further the one Class II (pOph-ISO 035) with weak Pa/3 
in absorption. 

The accretion luminosity of each Class II object is 
delved from the empirical correlation between L^cc and the 
luminosity of Pa/3 or Br7, derived by Natta et al. H2(JU4|I 
and Calvet et al. H2004|l . respectively (see also Muzerolle 
et al. nMH|l : 

log Lacc/i© = 1.36 \ogL{Pp)/LQ -h 4 (3) 

\ogLacclLQ = 0.9 (logi(Br^)/io + 4) - 0.7 (4) 

These relations have been calibrated using accretion 
luminosities derived by fitting the measured veiling (for 

http:/ /www. eso.org/instruments/isaac/tools/spectroscopic 
_standards.html 



T Tauri stars) and/or profiles with the predictions of 
magnetospheric accretion models; the objects used for the 
calibration cover the mass range from 3 Mq to brown 
dwarfs. The mass accretion rate is then computed from 
Lace {Mace = L acc R* / {G AU)) . Thc rcsults are given in 
Table C.l. 

The reliabilty of our procedure was verified by apply- 
ing it to a sample of well studied pre-main sequence stars, 
covering roughly the same range of masses, for which reli- 
able values of the stellar parameters (i.e., mass and radius) 
and of the accretion rate could be found in the literature. 
Using literature measurements of the Pa/3 intensity and 
of Lj,, we derived for each object mass and accretion rate 
as done for the Ophiuchus objects, and compared them to 
the "real" values. Details can be foimd in the Appendix 
B. 

We have applied a similar procedure to thc Class III 
objects; the results are shown in Table C.2. 

4. Results 

4.1. Equivalent Widths 

Emission in the near-IR hydrogen recombination lines has 
been detected in 45% of Class II sources, 46 of the 96 ob- 
served in Pa/3 and 1 out of 9 observed in Br7. In contrast, 
no Class III source shows emission in these hydrogen lines; 
8 Class III objects have Pa/3 or Br7 in absorption and for 
the others we do not see the lines. The measured equiva- 
lent widths are given in Table C.l and C.2. Fig. ^ shows 
the Pa/3 equivalent width as function of L^,. 

Six Class III objects have Pa/3 in absorption with 
equivalent widths >1 A, i.e., larger than one can expect in 
late- type stars (Wallace et al. I2000|l . They are likely ear- 
lier type stars, and this is certainly the case of p Oph-ISO 
180, which is classified A7 by Wilking et al. (2005) and of 
p Oph-ISO 113, earlier than F8 according to LR99. For 
these six stars, as already mentioned, the method used to 
estimate Aj and all the derived stellar parameters is not 
correct; therefore, we omit their stellar parameters from 
Table C.2. 

The comparison between the Class II and Class III 
samples clearly shows that emission in the near-IR hy- 
drogen lines, in contrast to that in optical lines such as 
Ha is restricted to objects with circumstellar disks, and 
can therefore be used as a reliable accretion indicator. 
However, one should keep in mind that the opposite is not 
necessarily true, as about 50% objects with disks have no 
detected emission. 

The fraction of Class II objects with detected Pa/3 
emission varies from 56% for L^>1 L0to 42% for 0.03 < 

< 1 Lq. Very low luminosity objects (7 objects with 
L^<0.03 L0) have a marginally higher detection rate 
(~57%), due in part to the sensitivity limit of our mea- 
surements, which is higher for lower luminosity objects 
(see Sec. 3.1), but also due to the incompleteness of the 
BKAOl survey for very low luminosity sources, which are 



Natta et al.: Accretion in p-Oph 



5 



detected only when they have a large mid-IR excess, very 
likely indicative of higher accretion rates. 
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Fig. 1. Equivalent width of the Pa/? emission line for all 
objects with known L^. Circles are Class II sources, filled 
circles are detections, empty circles with arrows upper 
limits. Squares are Class III sources. Two objects with 
Li,< 0.01 Lq and Pa/? undetected are not plotted. 



4.2. Accretion Luminosity 

Fig. 13 shows the accretion luminosity of Class II objects 
computed from the IR line luminosity as a function of L^,. 

For any given L^,, there is a large range of measured 
Lace (about 50), which does not seem to vary significantly 
with L^,; because of our sensitivity limit, this is probably 
just a lower limit to the actual range of Lace- One can 
also see that for the majority of objects Locc/L^<0.1, but 
there is a significant fraction of cases with Lacc~L^. 
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Fig. 2. Accretion luminosity from the IR lines as func- 
tion of for Class II objects. Dots show Lace mea- 
surements from Pa/? (filled: detections, empty: upper lim- 
its); diamonds measurements from Br7 (filled: detections, 
empty: upper limits); The dotted lines show the locus of 
Lacc/L*=0.01, 0.1 and 1, as labelled. Crosses are objects 
in Taurus (see text for references). 

5. Discussion 

The results summarized in Fig. |3| describe the accretion 
properties of the largest sample of Class II stars in any 
single star-forming region studied so far. The sample con- 
tains more than hundred objects with evidence of disks, 
and is complete in the mass interval from ~0.03 to about 3 
Mq. The corresponding accretion rates vary from ~ 10~^^ 
to ^ 10^® Mo/yr, with a strong dependence of Mace on 
(Maccoc M^-*=^°-2). For any M*, there is a large dis- 
persion of values of Mace, of two orders of magnitude at 
least, which does not seem to change with M*. Note that 
the real spread is likely bigger, because of the many upper 
limits in our survey. 



4.3. Mass accretion rate 

Fig. 121 shows the mass accretion rate Mace of Class II 
sources as function of M^,. There is a clear trend of in- 
creasing A'lacc with increasing M^,. Not including upper 
hmits, we find using ASURV (Feigelson and Nelson[l985) 
MaccO^ Mj,^-®^°-^; the slope does not change if we include 
the upper limits in the analysis. 

Superimposed on this trend, there is a large spread 
of Mace for any value of M^, of two orders of magnitude 
at least. Within statistical fluctuations, the objects are 
distributed quite uniformly in this range. 



5.1. Ophiuchus and Taurus 

One of the aims of our study was to compare the accretion 
properties in Ophiuchus with those of objects in Taurus. 
The Taurus results are shown by crosses in Fig. |21 and 
121 The accretion luminosity and mass accretion rate have 
been derived from the UV and optical veiling and/or by 
fitting with magnetospheric accretion models the Ha pro- 
file. This second method is the only possible one for very 
low mass objects and BDs, since veiling cannot be detected 
below a limiting value Macc^lO"^" Mq /yr. The results are 
from GuUbring et al. (p98j) . MuzeroUe et al. |(DM1IM71 
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Log M,/Mq 

Fig. 3. Mass accretion rate derived from the IR lines as 
function of M^,. Symbols as in Fig. [3 The dashed line 
shows the relation AfaccOcM^,^-®, derived from a statistical 
analysis using ASURV; the dotted line plots, for compar- 
ison, the relation MqccOcM^,. 

White fc Gh ez (IMll . White & Basri ll^Hmil . and 

Calvet et al. H2()()4|l : note that, for homogeneity, we have 
re-determined using DM98 tracks for all objects. 

The methods used to derive L^cc and Mace in the two 
regions are therefore different, since in Ophiuchus Mace is 
derived from the luminosity of the hydrogen recombina- 
tion lines. However, the relations (eq. 3 and 4) we used 
have been "calibrated" mostly using Taurus objects (see, 
e.g., MuzeroUe et al. Il998i Calvet et al. i2004> Natta et 
al. I2()()4|l . so that we do not expect any systematic dif- 
ference in the Ophiuchus- Taurus comparison due to the 
different methods. 

The two figures show that the accretion properties of 
the two star forming regions are very similar. MuzeroUe 
et al. ij^OSj) derive MaccOc M^,^-^ for their sample (mostly 
in Taurus, with additional brown dwarfs from other star- 
forming regions), neglecting upper limits. Within the 
errors, this relation is identical to what we obtain in 
Ophiuchus. If we concentrate on Fig. O we can see that 
not only the slope of the relation of Mace with M^,, but 
also the range of values is very similar. In particular, the 
two samples have similar values of the maximum Mace for 
any given M*, and similar spread of Mace values, at least 
for M^>0.06 - 0.08 M©. 

For lower M^,, most Taurus BDs have very low accre- 
tion rates, 1-2 orders of magnitude lower than similar ob- 
jects in Ophiuchus. As already discussed, the fact that we 
do not find these very low accretors in Ophiuchus most 



likely reflects the incompleteness of the BAKOl sample at 
very low masses, and selects objects with comparatively 
strong mid-IR fluxes. Natta et al. H2U(J2|I showed that the 
BAKOl sample of brown dwarfs has relatively large lu- 
minosity, and is probably very young. As discussed in 
Sec. 4.1, a fraction larger than for more luminous objects 
has detected Pa/3. All this indicates that there may be low 
Mace BDs which are missing from the Ophiuchus sample. 
It is, in a way, more surprising that very few, if any, of 
the brown dwarfs in Taurus have high Mace, while higher 
mass objects in the two regions have very similar accre- 
tion properties. It is possible that this difference between 
the two regions at the very low end of the M* distribu- 
tion contains important information, that needs further 
investigation. This is, however, beyond the scope of this 
paper. 

5.2. Variability 

All pre-main sequence stars are variable objects, and, in 
particular, all the accretion indicators in TTS and BDs 
show large variability. 

Variability does not affect the correlation of Mace with 
Mi, , as the p Oph sample is sufficiently large that individ- 
ual fluctuations cannot change it. It may be more impor- 
tant when we consider the spread of Mace values for any 
given M^,. Recently, Scholz & Jayawardhana H2()()5|l have 
studied the variability of accretion indicators (mostly H^) 
for six young brown dwarfs; they claim that the accretion 
rate in some of their objects varies by at least one order of 
magnitude, and that this variability may account for the 
large spread in the Mace- M^, correlation. 

We have estimated the magnitude of the spread in 
Mace for individual objects by looking at the results of 
Gatti et al. (2006), who have recently obtained J-band 
spectra of a small (14 objects) subset of our Ophiuchus 
sample. The Gatti et al. sample includes both TTS and 
BDs, observed one to two years later than the spectra dis- 
cussed in this paper. The two data sets show variations 
in the Pa/3 equivalent width of a factor of two at most 
(in both directions), with only one exception, where the 
Pa/3 equivalent width has increased by a factor of three 
over the time interval between the two sets of observa- 
tions. For the same objects, we have also looked in the 
literature for variations of the broad-band J magnitude, 
used to compute the line flux (Sec. 3. 2). The variation of 
Mace, computed taking the maximum variations in the J 
magnitude and in the Pa/3 equivalent width, is of a factor 
~ 4. This is much smaller than the dispersion of points 
in Fig. 121 and would not change signiflcantly any of our 
conclusions. 

A detailed analysis of the variability of the IR emission 
lines and continuum, in analogy to what has been done for 
Hq (e.g., Johns-KruU & Basri I19ll7| . is certainly needed. 
However, from the results obtained so far, it seems unlikely 
that the dispersion of Mace values can be accounted for by 
variability alone, and that, if averaged over a sufficiently 
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long period of time, one would find that all the p Oph 
stars of a given mass accrete at the same rate. 

5.3. Viscous disks 

The Mace dependence on is difficult to understand in 
terms of disk physics, as discussed, e.g., by MuzeroUe et 
al. ij^UnSjl . Natta et al. PMji . Calvet et al. (;2004). In 
a standard steady accretion disk model, Mace is propor- 
tional to the disk mass divided by the time scale for vis- 
cous evolution. In an a-disk (Shakura & Sunvaev 11978(1 . 
the viscosity depends on the ratio ^/cj., where il is the 
keplerian angular velocity and Cs the sound speed; then, 

— 1/2 

Maccoc Md X Mi, X Td, where Md and Td are disk mass 
and temperature, respectively. With the further assump- 
tions that Md oc Mi, (e.g., Natta et al. l2()()()|l . and that the 
disk heating is dominated by the stellar irradiation, this 

1/2 

gives, to zero order, MqccOC x Ti,. For PMS stars, 

the relation between Toff and is rather shallow (ap- 
proximately T* cx M° "* for M*>0.1 M0, and much flatter 
for lower masses; see, e.g., DM98) and we expect Mace to 
increase roughly as M2 , with 7^1. The relation will be 
even flatter if the contribution of the stellar radiation to 
the disk heating is negligible. 

It is possible that a (or, more generally, the efhciency 
of momentum transfer) depends, in turn, on M*. If viscos- 
ity is the result of magneto-rotational instabilities (MRI) 
(see, e.g., Balbus & Hawlev Fiyyifl . the disk gas should be 
sufficiently ionized. Muzerolle et al. (2003J suggest that 
the steep correlation of Mace with Mj, can be explained 
if the disk ionization is controlled by the X-ray radiation 
from the star, since the X-ray luminosity is not constant 
over the mass spectrum, but is observed to increase with 

X-ray observations of Ophiuchus have been recently 
carried out with Chandra and XMM satellites by Imanishi 
et al. I^Mt and Ozawa et al. (^UDF). Both studies de- 
tected a significant fraction of Class II sources (70 and 48 
% respectively) ; they found that the X~ray spectral prop- 
erties, as well as the relationship between Lx and Lboi of 
class II sources are similar to those of class III sources, 
but did not investigate the behaviour of X-ray luminos- 
ity with stellar mass. To our knowledge, the only study 
addressing the relation between mass and X-ray luminos- 
ity for young stars over a large range of luminosities and 
masses is in Orion. The COUP Chandra observations of 
Orion show that Lx scales approximately as M^^-^"^ "* in 
the interval 0.1-2 Mq (Preibisch et al. l20U5|l . However, it 
is not clear that this variation of Lx is sufficient to pro- 
duce the observed Macc-M*correlation, and more detailed 
MRI models, which include X-ray ionization, are required. 
If the X-ray emission of the central star is controlling ac- 
cretion, the large spread of Lx observed in the COUP 
data could also explain the large spread of Mace for any 
given M^. 

Viscous disk models predict that Mace decreases with 
time (e.g., Hartmannet al. ll998|l . Calvet et al. ( , 2000,1 esti- 



mate AlaccOc t~^'^, with a large uncertainty, from a sample 
of TTS in Taurus, Chamaeleon and Ophiuchus. Neither 
the similarity of accretion rates between Ophiuchus and 
Taurus nor the very large spread observed in both regions 
support age as a main factor in the determination of Mace- 
If the Calvet et al. H20()0|l rate is correct, the difference in 
age between Taurus and Ophiucus should give on average 
a difference in Mace of a factor >3, of which we have no 
evidence. In addition, the Ophiuchus Mace range of more 
than two orders of magnitude corresponds to an age range 
of at least a factor 20, much too large when compared to 
the HR location of the objects (see, e.g., LR99). 

The time evolution of viscous disks is influenced by 
the presence of close companions (see Calvet et al. I2()()()|l . 
Companions truncate the circumstellar disk at a radius 
which depends on the binary separation. As the disk 
evolves, more and more matter expands outside the trun- 
cation radius, with the effect of decreasing the disk mass 
and Mace- A sample of objects with the same initial value 
of Mace but companions at different distances will show 
with time an increasing spread of Mace values. 

This effect, however, is not seen in the Taurus TTS 
(White & Chez I^IHTTIl . where the accretion rate is similar 
for single and primary stars with companions as close as 
10 AU. At the age of Ophiuchus, only very close compan- 
ions have had time to reduce Mace by a significant factor 
(separation <30 AU or <0.2 arcsec for an age of 10^ years 
according to Calvet et al. 2000(1 . There have been a num- 
ber of multiplicity surveys of Ophiuchus, some capable of 
detecting very close binaries. Three Class II objects (i.e., 
objects with a mid-IR detected circumstellar disk) have 
companions closer than <0.25 arcsec (Barsony et al. l2UU5l 
Ratzka et al. 12005(1 : one has detected Pa/?, while in the 
other two cases the line has not been detected. The ob- 
servational evidence of a correlation between the accre- 
tion rate and the presence of very close companions is 
clearly inconclusive. At this stage, it cannot be quantita- 
tively confirmed nor dismissed, and should be investigated 
further. 

5.4. Initial conditions 

Although all the effects discussed so far can play a role and 
need further investigation, it is possible that differences in 
the initial conditions, i.e., in the physical properties of the 
molecular cores from which the star-fdisk system forms, 
determine the TTS disk properties, and in particular the 
behaviour of Alacc disussed in this paper. 

The self-similar viscous disk models of Hartmann et 
al. 1(1998(1 show that the accretion rate is proportional to 
the disk mass at t=0, i.e., when accretion onto the disk 
stops, and, in the early phases ot the evolution, to its t=0 
outer radius, which in turns depend on the core properties. 
Alexander & Armitage l(2()()f)(l have started exploring how 
this can introduce a MaccOC M*^ correlation at a later time. 

More realistic models that follow the formation and 
evolution of circumstellar disks (Hueso & Guillot 12005(1 
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illustrate clearly how different core properties (in particu- 
lar, different rotation velocities) can create a large spread 
of Mace for objects with the same and age. 

Models that compute the evolution of disks starting 
from the core infall phase over a large range of parame- 
ters are required, if we want to estimate the effect of the 
initial conditions on the relation of Mace with M^, and on 
its scatter. The observations presented in this paper, and 
the similar results for Taurus, provide an excellent test 
of such models. Note that the the fact that disk accre- 
tion properties in Taurus and Ophiuchus are very similar, 
while the two regions have large differences in their envi- 
ronment, should put strong constraints on these models, 
which will be interesting to explore fully. 

6. Summary and conclusions 

In this paper, we report the results of a near-IR spectro- 
scopic survey of a large sample of very young objects in the 
p Oph core. The sample includes all Class II objects, i.e., 
objects with evidence of circumstellar disks from mid-IR 
photometry (BKAOl). This sample covers the mass range 
between about 0.03 to 3 Mq; according to BKAOl, it is 
complete to a limiting magnitude of about 0.03 L©, or 0.05 
M0. We have also observed a significant fraction of Class 
III objects, i.e., with no mid-IR excess emission, covering 
a similar range of luminosities. 

In contrast to the Balmer lines, the near-IR hydrogen 
recombination lines are seen in emission only in Class II 
objects. Of all our Class III sample, none has detected Pa/3 
emission. This confirms our assumption (Natta et al.[2004) 
that the near-IR lines can provide an immediate indication 
of the accreting properties of young stars, even when only 
relatively low resolution spectra are available. 

We have derived the mass accretion rate Mace from the 
luminosity of the hydrogen recombination lines, mostly 
from Pa/3 but in few cases from Br7. In total, we obtain 
measurements of Mace for 45 Class II objects, and upper 
limits for 39. 

Our results show that Mace increases sharply with 
M*(cx; Mj-^^°-^). We also find a large range of values of 
Mace for any given value of M^, (a spread of roughly two 
order of magnitudes, independent of M*). As discussed in 
the text, this is likely a lower limit to the true dispersion. 

When compared to accretion measurements in Taurus 
(see Muzerolle et al. l2UU5l and references therein), we find 
that the two regions look very similar, at least for objects 
with M^,>0.1 Mq. For both Taurus and Ophiuchus, the 
dependence of Mace on M*, the upper envelope of the 
Mace distribution (i.e., the largest values of Mace that any 
object of a given mass seems able to sustain), and the 
range of Mace values for any given M^, are very similar. At 
lower mass, the accretion rates of the Ophiuchus objects 
are much larger than their Taurus analogs. 

The observed behaviour of Mace does not have an ob- 
vious explanation. The correlation of Mace with M^, may 
be due to a dependence of the disk physics on the proper- 
ties of the central star. Muzerolle et al. (,2003^) suggest as a 



cause the effect of the X-ray emission from the central star 
on the disk ionization and angular momentum transfer. It 
is also possible that the correlation reflects the properties 
of the pre-stellar cores, from which the star and disk form. 
Both possibilities need to be explored further. 

The large spread of values of Mace for any M^, may 
be a side-product of the same mechanisms that produce 
the correlation between these two quantities, as discussed 
in Sec. 5. In addition, other effects may play an important 
role, for example the dynamical action of close compan- 
ions, or the intrinsic variability of the accretion process. 
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Appendix A: Testing the assumption of coeval 
star formation 

The assumption of coeval star formation, albeit quite rea- 
sonable for a region like p Oph, introduces errors in our 
results. The same is true of the choice of any specific set of 
evolutionary tracks. However, it turns out that both kinds 
of errors are unimportant, when dealing with a large sam- 
ple of objects as in our case. 

Fig. lA. II shows the analog of Fig.O reproduced on the 
top-left panel, computed using the DM98 isochrone for 1 
My and the evolutionary tracks of Siess ( 20001 for 0.5 My 
and 1 My, respectively. Older tracks give slightly lower 
values of Mace, especially for more massive objects, while 
the range of M^ remains practically the same. Adopting 
different evolutionary tracks does not change the results. 
The main consequence of assuming coeval star formation 
is to reduce slightly the real spread of Mace for any given 
value of M^. 

Appendix B: Taurus: a test of the method 

A validation of the method used to compute the two quan- 
tities Mace and M^^and an estimate of the errors can be 
obtained by applying the same procedure to a sample of 
objects with known stellar parameters and accretion rates. 

The only sample for which this is possible is Taurus, 
which has been studied extensively over a large range of 
masses We have taken all the Taurus objects for which 
we could find in the literature reliable stellar parameters 
and accretion rates, measured from veiling and/or by fit- 
ting the observed Ha profiles with magnetospheric accre- 
tion models (Muzerolle et al. 1998, 2003, 2005; Calvet et 
al. I2004|l . For those with published Pa/3 fluxes or equiv- 
alent widths, we have followed the same procedure used 
for the p Oph sample. We have flrst computed Lace from 
L(Pa/3), and determined the stellar parameters M^/R^, 
and M* from L*, assuming coeval star formation at 1 My 
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Fig. A.l. Same as Fig.Ofor different ages and evolution- 
ary tracks. The top-left panel is for DM98 0.5My (as in 
Fig. 13, the bottom-left for DM98, 1 My, the top right is 
for Siess l|2()()()|l evolutionary tracks at 0.5 My, the bottom 
right at 1 My. In each Panel, the two dotted lines {Maccoc 
M*^) have been drawn to guide the eye in the comparison. 

and the DM98 evolutionary tracks. Mace is then derived 
from Lace and M^/R^. 

The results are summarized in the Fig. IB. II The top 
panel shows the complete sample of Taurus objects for 
which we could find measurements of Mace in the litera- 
ture. The squares are those for which also Pa/3 data exist; 
because none of the BDs in Taurus has a published J- 
band spectrum, we have added the BDs in Ophiuchus and 
Chamaeleon for which Natta et al. H2004|l have measured 
Mace from model fitting of the Ha profiles. The bottom 
panel shows the same plot when both M^, and Mace are 
derived from the observed L^, and Pa/? luminosity, as for 
the p Oph stars. 

The results indicate that our procedure does not intro- 
duce systematic trends in the results. The trend of Mace 
increasing sharply with M^, is reproduced in our method, 
and also the range of Maec for a given M* is similar, even 
if, as expected, the assumption of coeval star formation 
underestimates its spread slightly. 

Appendix C: Tables 

Tables are available in electronic form only. 
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Table C.l. Class II Objects 



(1) 


(2) 








(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 


(15) 


(16) 


(17) 


# 


Object 










Coordinates 


J 


H 


K 


Aj 


Lg U 


Lg Te// 


Lg M* 


EW (Pa,3) 


Inst. 


L(Pa^j) 


Lace 


Mace 


Other 




(ISO# 


) 








(J2000.0) 




(mag) 




(mag) 


(Lo) 


(K) 


(Mo) 


(A) 




(Lo) 


(Lo) 


(Mo/y) 


Names 


1 


pOph- 


[SO 


001 


16 25 36.74 


-24 15 42.40 


10.42 


9.04 


8.38 


1.8 


0.21 


3.59 


-0.26 


<-0.7 


S 


<-4.27 


<-1.81 


<-8.66 


IRS2 


2 


pOph- 


[SO 


002+ 


16 25 38.12 


-24 22 36.30 


12.84 


10.75 


9.54 


3.1 


-0.23 


3.53 


-0.52 


-2.8 


S 


-4.07 


-1.54 


-8.22 


B162538- 242238 


3 


pOph- 


[SO 


003 


16 25 39.58 


-24 26 34.90 


11.89 


10.05 


8.95 


2.5 


-0.10 


3.55 


-0.44 


-3.9 


s 


-3.81 


-1.18 


-7.92 


IRS3 


4 


pOph- 


[SO 


006 


16 25 56.16 


-24 20 48.20 


9.15 


8.14 


7.52 


0.5 


0.18 


3.59 


-0.28 


-19.0 


s 


-2.87 


0.10 


-6.74 


SR4/IRS12 


5 


pOph- 


[SO 


009 


16 26 1.37 


-24 25 20.40 


14.43 


12.44 


11.24 


2.8 


-1.07 


3.45 


-1.03 


<-0.3 


I 


<-5.74 


<-3.80 


<-10.24 


SKSl-4 


6 


pOph- 


[SO 


012 


16 26 4.58 


-24 17 51.50 


15.79 


13.41 


12.19 


4.2 


-1.05 


3.45 


-1.02 


<-0.2 


I 


<-5.96 


<-4.11 


<-10.55 


B162604-241753 


7 


pOph- 


[SO 


013 


16 26 7.04 


-24 27 24.20 


15.35 


12.38 


10.64 


5.2 


-0.38 


3.51 


-0.62 


-2.1 


s 


-4.34 


-1.90 


-8.53 


B162607-242725 


8 


pOph- 


[SO 


017 


16 


26 


10.33 


-24 20 54.80 


14.37 


10.85 


8.47 


5.8 


0.33 


3.61 


-0.19 


<-0.8 


s 


<-4.11 


<-1.58 


<-8.48 


GSS26 


9 


pOph- 


[SO 


019 


16 


26 


16.84 


-24 22 23.20 


11.03 


9.13 


8.20 


3.1 


0.54 


3.64 


-0.06 


<-0.9 


s 


<-3.87 


<-1.26 


<-8.23 


GSS29/EL18 


10 


pOph- 


[SO 


020 


16 


26 


17.06 


-24 20 21.60 


9.65 


8.61 


8.06 


0.8 


0.10 


3.58 


-0.33 


-1.2 


s 


-4.15 


-1.65 


-8.46 


DoAr24/GSS28 


11 


pOph- 


[SO 


023 


16 


26 


18.82 


-24 26 10.50 


14.84 


13.20 


12.14 


1.8 


-1.69 


3.42 


-1.40 


-1.8 


I 


-5.59 


-3.61 


-9.93 


SKS1-BDN04 


12 


pOph- 


[SO 


024 


16 


26 


18.87 


-24 28 19.70 


12.58 


9.93 


8.07 


3.7 


0.18 


3.59 


-0.28 


-8.9 


I 


-3.20 


-0.35 


-7.19 


VSSGl 


13 


pOph- 


[SO 


026 


16 


26 


20.97 


-24 8 51.90 


10.88 


9.87 


9.50 


1.0 


-0.33 


3.52 


-0.59 


-2.4 


I 


-4.23 


-1.76 


-8.41 


RBR15 


14 


pOph- 


[SO 


030 


16 


26 


21.53 


-24 26 1.00 


12.57 


11.52 


10.92 


0.7 


-1.22 


3.44 


-1.12 


-0.3 


I 


-5.95 


-4.09 


-10.50 


GY5 


15 


pOph- 


[SO 


032 


16 


26 


21.90 


-24 44 39.80 


12.34 


11.48 


10.86 


0.0 


-1.43 


3.43 


-1.24 


-0.4 


I 


-6.01 


-4.17 


-10.54 


GY3 


16 


pOph- 


[SO 


033 


16 


26 


22.27 


-24 24 7.10 


16.45 


15.09 


13.94 


0.6 


-2.95 


3.37 


-2.16 


<-0.7 


I 


<-7.15 


<-5.73 


<-11.81 


GYll 


17 


pOph- 


[SO 
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1 p; fin 

lo.oy 


1 o m 
Iz.U / 


n OK 

y.yo 


f; 1 
0. / 


n 1 f; 
U.IO 


Q p; Q 

O.OO 


n on 
-U.zVJ 


Q 7 

-o. / 


b 


Q ftn 
-O.OU 


n Qn 

-U.oy 


7 70 

- / . /z 


O VOftO 

U Y zOz 


/U 


pU pa- 


LbU 14z 


1 OI OI QQ 

10 z / z / .OO 


o /I Q 1 If; f;n 
-z4 ol 10. OU 


1 O Q 

Iz.oo 


1 n QQ 
lU.oo 


n QO 

y.oz 


Q n 
o.U 


n n7 
-U.U / 


Q p; r: 
O.OO 


n AO 
-U.4o 


^ n p: 


b 


^ A ai 
<-4.0 / 


^ O Q 

<-z.oo 


^ o 1 n 

<-y.iu 


1 rcjCOOf^ /OVOft7 

VbbUzo/U YzO/ 


/ 1 


rt/^ T-iU 
pUpn- 


rc/^ 1/1/1 
LbU 144 


1 a OI OQ /I 
10 z/ zo.4o 


OA Ol Ol nn 
-z4 z / zl.UU 


1 p; iA 
lo. i 4 


1 O Q1 

Iz.ol 


1 n 1 n 
lU.lU 


c; Q 
O.o 


n Ofi 
-U.zO 


Q p:q 

O.OO 


n K A 
-U.04 


^ 1 /t 
<-1.4 


T 
1 


^ A An 
<-4.4U 


^ -1 fiQ 

<-i.yo 


^ Q ftft 

<-o.00 


TT^Q/l P; /f~^"V07Q 

lKb4o/U Yz/o 


TO 


pUpn- 


C/^ 1 A1 

LbU 14 ( 


1 fi 07 Qn 1 Q 

10 z^ oU.lo 


O/l 07 AO An 

-z4 z ( 4o.4U 


1 p; QO 
lo.oz 


1 1 p;o 
11. oz 


n no 

y.uz 


fi fi 

0.0 


n 07 
U.z / 


Q fin 
O.OU 


n OO 
-U.zz 


^ n PC 
<-U.o 


b 


^ A Qft 

<-4.oO 


^ 1 HQ 

<-i.yo 


^ Q QI 

<-o.ol 


TT>0/17 /OV070 

lKb4/ /Ijt Yz/y 


/ O 


pUpn- 


LbU lol 


1 fJ 07 *3n Q A 

10 z/ oU.o4 


OA OA ccf; nn 
-z4 z4 oO.UU 


1 o in 
IZ. /U 


1 n ncc 

lu.yo 


1 n ni 
lU.U/ 


o f; 
Z.O 


n /in 
-U.4U 


Q PCI 
O.Ol 


n fio 
-U.Oz 


^ n PC 
<-U.O 


b 


<-4.y / 


^ o la 
<-z. /O 


^ n Qn 

<-y.oy 


OVOQ/1 
Vjr Y zo4 


TA 

( 4 


pUpil- 


■Q/^ 1 t^A 

LbU l04 


1 f; 07 QO Q f: 
10 z/ oz. OO 


O/I QO OA Qn 

-z4 oz o4.oU 


1 fi 1 n 

10. ly 


1 O 1 A 

LZ. { 4 


1 n nfi 

lu.yo 


fi n 

O.y 


n ni 
-U.Ul 


Q Plfi 
O.OO 


n QO 
-U.oy 


^ 1 n 
<-i.U 


b 


<-4.oz 


^ 1 QQ 

<-l.oo 


^ Q ftPC 

<-o.0o 


O V001 

U Y zyi 


/O 


pUpn- 


LbU loo 


1 a OI QQ 1 1 

Id Z{ OO. 11 


OA /II ir^Qn 
-z4 41 lO.oU 


1 1 QO 

11. oz 


n 1 Q 

y.io 


7 O 1 


Q o 
O.Z 


n /I n 

u.4y 


Q fiQ 
O.DO 


n nn 

-u.uy 


A n 
-4.U 


c 
b 


Q Oft 

-o.zO 


n AO 
-U.4o 


7 Qn 

- /.oy 


r^\zono 

Kj Y ZyZ 


( 


pvjpn- 


CO 1 fin 

LoKJ IDU 


1 fi 01 /IO 
ID Z / ( .4Z 


0/1 1 7 c,/i on 
-Z4 1 i 04. yu 


1 /I 1 

14.10 


1 7fi 
IZ. / 


1 1 op; 
1 i.yo 


1.0 


1 P.7 
-1.0/ 


/IQ 
0.40 


1 

-1.00 


Q Q 
-0.0 


T 
1 


p; 00 

-O.ZZ 


Q HQ 

-o.uy 


Q AA 
-y.44 


R1 ft07Q7 0/1 1 7P;ft 
I310Z / /-Z41 / 00 


i ( 


pUpn- 


OO 1 fiQ 

LbU ioo 


1 a. OI QQ QO 

10 z< oo.oz 


O/I Qfi c:q fin 
-z4 oO OO.OU 


1 1 QQ 
11. OO 


n AO 

y.4o 


Q 07 
O.Z i 


O 1 
Z. { 


n OQ 
U.zo 


Q cn 

o.oy 


n OK 
-U.zo 


o n 
-z.U 


b 


Q on 
-O.oU 


1 1ft 
-1.10 


Q no 
-o.Uz 


TT>0/in /OVQnQ 

lKb4y / U Y oUo 


1Q 

to 


pUpn- 


LbU io4 


1 f; o7 QQ f;Q 
10 z^ oo.Oo 


OA QQ Qn on 
-z4 OO oy.zU 


1 Q 07 

lo.z ( 


1 1 nQ 

11. yo 


1 1 no 
11. Uo 


1 o 
i.z 


1 Q1 

-l.oi 


OAA 

0.44 


117 
-1.1 / 


n Q 
-U.O 


T 
1 


PC fin 
-O.OU 


Q fil 
-O.Ol 


1 n ni 
-lU.Ul 


f^VQi n 
Y olU 


Id 

( y 


pUpn- 


CO 1 p: 
LbU IDO 


1 f? 07 QQ n/1 

10 z/ oo.y4 


OA An on in 
-Z4 4U zU. / U 


1 fi ^ A 

10.04 


1 Q ni 

lo.yi 


1 o on 
Iz.zy 


A O 

4.Z 


1 Q p: 
-l.oo 


O A A 

0.44 


1 on 
-l.zU 


n ft 

-y.o 


T 
1 


A Plft 

-4.00 


o on 
-z.zU 


Q PCO 
-O.OO 


O VQI O 

U Y olz 


ou 


pwpii- 


qo 1 ftft 

LOW lUU 


1 ft 97 QQ 4Q 
ID z / oy.^o 


OA '^Q c,n 
-Z'i oy ±o.ou 


lU. / o 


Q 01 
y .z ± 


0.4:U 


9 9 
z.z 


n 94 


o.oy 


n 9^ 
-u.zo 


Q Q 
-o.o 


T 
1 


-O.OO 


n f^7 

-U.O / 


7 7Q 
- ( . / o 


V 1 1 014: 


81 


pOph- 


[SO 168 


16 27 40.29 


-24 22 4.00 


8.44 


7.67 


7.21 


0.0 


0.28 


3.60 


-0.22 


-1.3 


S 


-3.96 


-1.38 


-8.26 


SR9/GY319/IRS52 


82 


pOph- 


[SO 170 


16 27 41.61 


-24 46 44.70 


17.20 


15.33 


13.55 


1.0 


-3.07 


3.36 


-2.23 


<-3.0 


I 


<-6.63 


<-5.01 


<-11.07 


B162741-244645 


83 


pOph- 


[SO 171 


16 27 41.75 


-24 43 36.10 


0.00 


14.88 


12.29 


8.5 


-0.50 


3.50 


-0.68 


<-4.0 


I 








GY323 


84 


pOph- 


[SO 172 


16 27 42.70 


-24 38 50.60 


13.24 


11.44 


10.54 


2.8 


-0.56 


3.49 


-0.72 


<-1.0 


s 


<-4.82 


<-2.56 


<-9.12 


GY326 


85 


pOph- 


[SO 175 


16 27 45.79 


-24 44 53.90 


17.38 


14.54 


12.46 


4.0 


-1.82 


3.42 


-1.48 


-15.7 


I 


-4.77 


-2.49 


-8.79 


GY344 



Table C.l. continued. 



oo 


pUpn-ibU i (u 


1 a OT Ad on 


o/i 01 /II on 
-z4 ol 41. zU 


1 Q OO 


1 O Ol 

Iz.zl 


1 1 OO 
11. OZ 


O O 
Z.Z 


1 11 
-1.11 


OAK 

o.4o 


1 nc 
-l.Uo 


^ n c 
<-U.O 


1 


^ C £!0 


^ O £i A 

<-o.d4 


^ 1 n no 
<-lU.Uo 


/^"vocn 
(orYooU 


O ( 


pUpii-ioU Vi i 


1 07 /1 7 no 


0/1 in 

-z4 40 oo.lU 


lo. / 


1 O Ol 

Iz.ol 


11 1 O 

11. lo 


c: o 
O.z 


n kp: 
-U.OO 


o /in 
o.4y 


n 70 
-U. / z 


^ n /I 


T 
1 


^ ^ Ol 

<-o.zl 


^ n nn 


<-y.DD 


Kj Y OOZ 


o o 


pUpn-lbU i (o 


1 i5 07 /I n 70 


0/1 oc oo nn 
-z4 ZD zz. UU 


1 O 70 

Iz. 


11 1 o 

11. Iz 


1 n 1 i5 
lU.lD 


o o 

Z.Z 


n i5 c 
-U.DO 


O AO 


n 77 
-U. / / 


o n 
-2.U 


C! 

b 


A cir\ 
-4.dU 


O O K 
-Z.ZO 


O 70 

-o. to 


i^\/071 


89 


„/^„i TC/^ IOC 

pUpn-lbU 185 


"1 £! 07 C C O C 

Id z7 55.25 


o /I oo on /?n 

-z4 a6 o9.dU 


13.04 


11.59 


10.79 


1 7 

1.7 


-0.98 


0.45 


-0.97 


<-0.2 


T 
1 


<-5.90 


<-4.03 


<- 10.48 


/^■\/"on7 


on 


pUpil-iow lo/ 


1 R 07 f^r^ 
Id z/ 00. oo 


0/1 ofi 1 1 on 
-z4 ZD i / .yu 


iU. 14 




E5.yu 


U.z 


-u.oy 


o.Ol 


-U.Dz 


-O.U 


Q 

O 


Q Q1 

-o.yi 


-l.OZ 


- / .yo 


blxlU/ Lt Y 4UU 


91 


pOph-ISO 190 


16 28 3.56 


-24 34 38.60 


0.00 


15.24 


13.20 


6.1 


-1.37 


3.44 


-1.21 


<-5.0 


I 








GY450 


no 
yz 


^r^^u TC/^ 1 no 
pUpil-ibU lyo 


iO 28 iz. / Z 


O/i 11 occc?n 
-z4 11 OO.DU 


1 O 1 

lo.Dl 


1 o no 
Iz.Uz 


1 1 nn 

11. uy 


1 n 

i.y 


1 11 
-1.11 


O /I K 

o.4o 


1 nc: 
-l.Uo 


1 O 

-l.o 


T 
1 


K n7 
-O.U/ 


o on 

-z.oy 


n oo 

-y.oz 


"Dlt^OOlO 0/I1100 

J3lDZolz-z411oo 




pwpii low 


16 28 13.79 


-24 32 49.40 




±\J.OiJ 


10 10 

J-U. -LU 


1 s 

i-.O 


-0 fi4 




-0 77 


^-0 ^ 


O 


^ ^ 1Q 


Ofi 


^-Q fiO 


JJ±UZO-LO Z^OZ^c/ 


94 


pOph-ISO 195 


16 28 16.73 


-24 5 14.30 


10.98 


9.57 


8.86 


1.8 


-0.04 


3.56 


-0.41 


-0.8 


s 


-AAA 


-2.03 


-8.79 


ISO 1628168-240519 


95 


/lOnh-TSO 196 


16 28 16.51 


-24 36 58.00 


11.31 


10.08 


9.32 


0.9 


-0.55 


3.49 


-0.72 


-3.2 


s 


-4.31 


-1.86 


-8.43 


WSRfifl /Rl 6981 6-943657 


96 


rtOnh-TSO 199 


16 28 45.60 


-24 28 19.00 


9.21 


8.41 


8.00 


0.2 


0.02 


3.57 


-0.38 


-1.7 


s 


-4.06 


-1.53 


-8.31 


SR13 


41 




Coordinates 


T 
J 


TT 

H 


IV 


A, 7 


bg L* 




Ijg iVl* 




inst. 




T 

-L'acc 


Mace 


vjtiier 






(J2000.0) 




(mag) 




(mag) 






1^0 j 










/'l\/r /lA 


Names 


y ( 


pUptl-lbU uoy 


16 26 31.04 


-24 31 5.20 


14. 9d 


1 O QO 

Iz.oz 


1 n Qf; 


4.D 


-U.4y 


O.OU 


-U.DO 


<-U.t) 


b 


<-4.UD 


^1 K A 

<-1.04 




WL/ / /Lr Y yo 


yo 


pL/pn-ioiJ U<U 


16 27 41.61 


-24 46 44.70 


1 7 on 
1 / .zU 


10. OO 


lo.oo 


l.U 


-o.U / 


o.oD 


O OQ 
-Z.ZO 


^ 1 o 

<-l.z 


O 


<-4.UD 


<-o.U0 


^ O 1 

<-y.iz 


VV LiZ/\ji 1 Izo 


99 


p(Jpn-lB(J U(0 


16 26 51.97 


-24 30 39.50 


U.UU 


1 a CO 


io.4D 


1 n y1 

1U.4 


-U.D9 


o.4( 


-U.8U 


<-o.9 


o 

b 


<-4.UD 


<-U.4U 


<-D.9z 


KjY 144 


100 


pOph-ISO 076 


16 26 53.47 


-24 32 36.20 


0.00 


16.24 


13.12 


10.7 


-0.48 


3.50 


-0.67 


<-3.0 


S 


<-4.06 


<-0.32 


<-6.92 


GY146 


101 


pOph-ISO 085 


16 26 58.28 


-24 37 41.00 


0.00 


0.00 


14.41 


-1.0 


-0.48 


3.50 


-0.67 


<-8.0 


s 








CRBR51 


102 


pOph-ISO 093 


16 26 31.04 


-24 31 5.20 


0.00 


15.65 


12.56 


10.6 


-0.24 


3.53 


-0.53 


<-0.7 


s 


<-4.06 


<-0.73 


<-7.41 


GY188 


103 


pOph-ISO 108 


16 27 9.43 


-24 37 18.80 


16.79 


11.05 


7.14 


10.5 


1.44 


3.76 


0.48 


-1.3 


s 


-4.06 


1.46 


-5.84 


EL29/GY214 


104 


pOph-ISO 139 


16 27 26.29 


-24 42 46.10 


0.00 


15.18 


12.66 


8.2 


-0.73 


3.46 


-0.82 


<-2.8 


s 


<-4.06 


<-0.77 


<-7.27 


GY260 


105 


pOph-ISO 161 


16 27 37.25 


-24 42 38.00 


0.00 


14.52 


11.46 


10.5 


0.27 


3.60 


-0.22 


<-0.8 


s 


<-4.06 


<-0.30 


<-7.17 


GY301 



Table C.2. Class III Objects 



(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 


(15) 


(16) 


(17) 


# 


Object 


Coordinates 


J 


H 


K 


Aj 


Lg L* 


Lg Te// 


Lg M« 


EW (Pa/s) 


Inst. 


L(Pa/3) 


^acc 


Mace 


Other 




(ISO#) 


(J2000.0) 




(mag) 




(mag) 


(Lg) 


(K) 


(M©) 


(A) 




(Lo) 


(Lo) 


(Mo/y) 


Names 


1 


pOph-ISO 005 


16 25 50.53 


-24 39 14.50 


9.98 


8.82 


8.33 


1.4 


0.21 


3.59 


-0.26 


<-0.5 


S 


< -4.42 


< -2.01 


<-8.86 


IRS 10 


2 


pOph-ISO Oil 


16 26 3.29 


-24 17 46.50 


10.67 


9.58 


9.12 


1.2 


-0.18 


3.54 


-0.49 


<-0.8 


I 


< -4.57 


< -2.22 


<-8.92 


VSSG19 


3 


pOph-ISO 014 


16 26 7.64 


-24 27 41.40 


14.68 


11.85 


10.41 


5.4 


-0.03 


3.56 


-0.40 


<-0.5 


s 


< -4.64 


< -2.31 


<-9.07 


B162607-242742 


4 


pOph-ISO 016 


16 26 9.31 


-24 34 12.10 


7.74 


6.95 


6.50 










-h5.7 


s 








SR3 


5 


pOph-ISO 018 


16 26 15.81 


-24 19 22.10 


14.03 


11.40 


10.03 


4.8 


-0.02 


3.56 


-0.40 


<-0.5 


s 


< -4.63 


< -2.29 


<-9.06 


SKSl-7 


6 


pOph-ISO 028 


16 26 21.02 


-24 15 41.50 


12.78 


10.38 


9.27 


4.5 


0.41 


3.62 


-0.14 


<-0.6 


s 


< -4.16 


< -1.66 


<-8.58 


B162621-241544 


7 


pOph-ISO 044 


16 26 28.48 


-24 15 41.20 


15.31 


12.37 


10.78 


5.4 


-0.29 


3.52 


-0.56 


<-0.7 


s 


< -4.73 


< -2.43 


<-9.09 


B162628-241543 


8 


pOph-ISO 064 


16 26 43.76 


-24 16 33.30 


12.98 


10.76 


9.60 


3.8 


0.00 


3.56 


-0.38 


<-0.7 


s 


< -4.46 


< -2.07 


<-8.84 


VSSGll 


9 


pOph-ISO 066 


16 26 44.30 


-24 43 18.00 


10.99 


10.02 


9.57 


0.7 


-0.51 


3.49 


-0.69 


<-0.8 


s 


< -4.87 


< -2.62 


<-9.20 


GY112 


10 


pOph-ISO 069 


16 26 47.05 


-24 44 29.90 


12.33 


11.12 


10.56 


1.3 


-0.83 


3.46 


-0.89 


<-0.4 


s 


< -5.47 


< -3.44 


<-9.92 


GY122 


11 


pOph-ISO 073 


16 26 49.23 


-24 20 2.90 


12.20 


9.85 


8.69 


4.2 


0.53 


3.63 


-0.07 


<-0.4 


s 


< -4.23 


< -1.75 


<-8.72 


VSSG3 


12 


pOph-ISO 074 


16 26 51.12 


-24 20 50.50 


13.81 


11.48 


10.21 










+1.5 


s 








IRS20/GY143 


13 


pOph-ISO 082 


16 26 56.92 


-24 28 37.10 


17.46 


14.75 


12.81 


3.8 


-1.94 


3.41 


-1.55 


<-2.0 


I 


< -5.77 


< -3.85 


< -10.12 


GY163 


14 


pOph-ISO 091 


16 27 1.62 


-24 21 37.00 


14.25 


11.07 


9.39 


6.1 


0.49 


3.63 


-0.09 


+0.9 


s 








VSSG8/GY181 


15 


pOph-ISO 113 


16 27 11.68 


-24 23 42.00 


14.21 


11.62 


10.11 


~ 








-1-4.0 


s 








IRS32/GY228 


16 


pOph-ISO 114 


16 27 11.71 


-24 38 32.10 


0.00 


15.06 


11.06 


16.2 


1.71 


3.79 


0.64 


<-8.0 


I 








WL19/GY227 


17 


pOph-ISO 135 


16 27 22.91 


-24 17 57.40 


13.33 


10.76 


9.45 










-hl.3 


s 








WSSG22 


18 


pOph-ISO 148 


16 27 31.06 


-24 34 3.20 


13.43 


11.36 


10.39 


3.6 


-0.27 


3.53 


-0.55 


<-0.5 


s 


< -4.86 


< -2.61 


<-9.28 


GY283 


19 


pOph-ISO 152^ 


16 27 32.68 


-24 33 23.90 


16.15 


12.74 


10.90 


6.6 


-0.10 


3.55 


-0.45 


<-2.0 


s 


< -4.10 


< -1.58 


<-8.31 


GY289 


20 


pOph-ISO 156 


16 27 35.26 


-24 38 33.40 


11.28 


10.23 


9.67 


0.8 


-0.62 


3.48 


-0.76 


<-0.4 


s 


< -5.27 


< -3.17 


<-9.71 


GY295 


21 


pOph-ISO 158 


16 27 36.52 


-24 28 33.30 


11.99 


11.38 


11.15 


0.0 


-1.29 


3.44 


-1.16 


<-1.0 


s 


< -5.48 


< -3.46 


<-9.86 


GY297 


22 


pOph-ISO 169a 


16 27 41.49 


-24 35 37.70 


14.30 


11.83 


10.56 


4.4 


-0.30 


3.52 


-0.57 


-1-0.7 


s 








GY322 


23 


pOph-ISO 169b 


16 27 41.64 


-24 35 41.10 


14.74 


12.48 


11.26 


3.8 


-0.78 


3.46 


-0.85 


<-0.9 


s 


<-5.06 


<-2.89 


<-9.38 


GY322 


24 


pOph-ISO 180 


16 27 49.87 


-24 25 40.20 


9.44 


8.12 


7.30 


~ 


~ 






-1-6.3 


s 








VSSG14/GY372 


25 


pOph-ISO 181 


16 27 50.51 


-24 39 3.10 


14.47 


12.74 


11.85 










+5.5 


s 








GY373 


26 


pOph-ISO 183 


16 27 51.92 


-24 46 29.60 


14.05 


11.61 


10.37 


4.3 


-0.22 


3.53 


-0.52 


<-0.7 


s 


< -4.66 


< -2.34 


<-9.03 


GY377 


27 


pOph-ISO 186 


16 27 55.65 


-24 44 50.90 


12.34 


11.15 


10.47 


1.0 


-0.98 


3.45 


-0.97 


<-1.5 


s 


< -5.03 


< -2.84 


<-9.29 


GY398 


28 


pUpn-ibO 189 


16 27 57.87 


-24 36 2.20 


15.37 


12.98 


11.83 


4.4 


-0.79 


3.46 


-0.86 


<-1.5 


c 
b 


< -4.86 


< -2.61 


<-9.10 


Cjr Y41z 


OQ 


/iOt-vV> TQO 1 QO 

pyjpii-ioyj lyz 


16 28 5.78 


-24 33 55.00 


ID. / 


1 A nn 


IZ.OO 


0.1 


1 ns 

-l.Uo 








T 

i 


^ A QQ 


^ 9 7Q 




V/1 70 


oU 


pUpn-lbU 19/ 


16 28 21.71 


-24 42 47.10 


1 C QQ 
10. OO 


14. Uo 


lz.44 


4.9 


1 1 Q 

-l.io 


o.4o 


-i.lU 


^ 1 o 

<-l.U 


T 

1 


< -o.oy 


< -o.oz 


^ O TA 

<-y. / 4 


r$lDzozl-z44z4D 


oi 




16 28 32.66 


-24 22 44.90 


O. lO 


/ .40 


D.?5D 


l.o 


U. / D 


O.D / 


U.U/ 


<-U.4 


c 


^ A flO 

< -^.{JZ 


^ 1 AT 
< -1.4/ 


^ Q CO 

<-O.J3Z 


or? on 


# 


Ubject 


Coordinates 


T 

J 


TT 
11 


K 


A 

A, J 


T ^ T 

Lg L, 


Lg T-eff 


T ^ 1\ T 

Lg IVl* 


tiW/ (DTf) 


Inst. 


L[DT-y) 


T 

J-'acc 


-'-''-'-ace 


Utner 




(ISO#) 


(J2000.0) 




(mag) 




(mag) 


(Lo) 


(K) 


(Mo) 


(A) 




(Lo) 


(Lo) 


(Mo/y) 


Names 


32 


pOph-ISO 047 


16 26 31.04 


-24 31 5.20 


14.96 


12.32 


10.86 


4.6 


-0.49 


3.50 


-0.68 


<-1.0 


s 


< -4.71 


< -1.34 


<-7.93 


IRS14/GY54 


33 


pOph-ISO 113 


16 27 11.68 


-24 23 42.00 


14.21 


11.62 


10.11 










+5.6 


s 








IRS32/GY228 


34 


pOph-ISO 157 


16 27 35.67 


-24 45 32.62 


12.71 


11.47 


10.88 


1.4 


-0.95 


3.46 


-0.96 


<-1.0 


s 


< -5.54 


< -2.08 


<-8.54 


GY296 


35 


pOph-ISO 179 


16 27 49.97 


-24 44 17.00 


13.85 


11.97 


10.94 


2.8 


-0.84 


3.46 


-0.89 


<-0.7 


s 


< -5.37 


< -1.94 


<-8.42 


GY370 



Caption of Tables EH and 
Column 1: running number; 

Column 2: ISOCAM number from BKAOl; a + sign marks objects with a companion not resolved in the 2MASS 
photometry; 

Column 3 and 4: J2000 coordinates 
Column 5,6,7: 2MASS J, H, K photometry; 
Column 8: J-band extinction; 

Column 9,10,11: stellar luminosity, effective temperature and mass, determined as described in the text; 

Column 12: line equivalent width: negative values for emission lines; 

Column 13: instrument used in the observations: S=SOFI/NTT, I=ISAAC/UT1; 

Column 14: line luminosity; 

Column 15,16: accretion luminosity and mass accretion rate; 
Column 17: other names. 



Table C.3. Companions not resolved by 2MASS 



(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


# 


Object 


Coordinates 


Separation 


EW 


Line 


Other 




(ISO#) 


(J2000.0) 


(arcsec) 


(A) 




Names 


1 


pOph-ISO 002b 


16 25 38.12 


-24 22 36.30 


1.8 


<-1.0 


Pa,3 


B162538-242238 


2 


pOph-ISO 068b 


16 26 46.43 


-24 12 0.10 


3.6 


<-0.3 


Pa/3 


VSS27 


3 


pOph-ISO 070b 


16 27 41.61 


-24 46 44.70 


4 


+2.5 


Br-y 


WL2/GY128 


4 


pOph-ISO 072b 


16 26 48.98 


-24 38 25.20 


3.6 


<-1.0 


Pa^ 


WL18/GY129 


5 


pOph-ISO 152b 


16 27 32.68 


-24 33 23.90 


3 


<-4.0 


Pa/3 


GY289 



